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1. Introduction 

The lattice simulation with the overlap fermions [|l|] provides a theoretically ideal setup to 
study low energy hadron physics. With the exact chiral symmetry, the continuum chiral perturba- 
tion theory (ChPT) can be applied without modification in sharpe contrast to other fermion formu- 
lations that violate either chiral or flavor symmetry. For this advantage, we performed numerical 
simulation using overlap fermion action with two dynamical flavors In particular, we carried 
out the calculation of light meson spectrum and studied the consistency between QCD and chiral 
perturbation theory. In this article, we present an extension of this study to the Nf = 2 + 1 QCD 
similar parameters. On a 16^ x 48 lattice, we generate 2,500 trajectories [Q, at ten different com- 
binations of up/down and strange sea quark masses, i.e. five rriu/s and two m/s. As in the Nf = 2 
case, topological charge of the gauge configuration is fixed to zero throughout Monte Carlo up- 
dates. Lattice spacing is determined by the Sommer scale as a = 0.1075(7) fm or a^' = 1.833(12) 
GeV from an input ro = 0.49 fm. 

In section |2[ we present basic part of the calculation of pseudo-meson masses and decay con- 
stants. Since these results depend on strange quark mass, the chiral extrapolation should be per- 
formed with a fit ansatz valid beyond the scale of kaon mass. To discuss this issue, we review the 
test of ChPT performed for the Nf = 2 case in section ||. Based on this test, in section ||, we present 
the extrapolation of the Nf = 2 + 1 data by using the two-loop ChPT formulae. 



2. Spectrum calculation (Nf = 2-1-1) 

For the Nf = 2+1 case, we calculate 80 pairs of the lowest-lying eigenmodes on each gauge 
configuration and store them on the disks. These eigenmodes are used to construct the low-mode 
contribution to the quark propagators. The higher-mode contribution is obtained by conventional 
CG calculation with significantly smaller amount of machine time than the full CG calculation. 
Those eigenmodes are also used to replace the lower-mode contribution in the meson correlation 
functions by that averaged over the source location (low-mode averaging) ^. The noise of 
correlation function is decreased as shown in Figure [I], where effective pion masses are compared 
between the data with and without the low-mode averaging for the lightest four quark masses. 

In order to obtain renormalization factor of quark mass Z,„, we calculate scalar and pseudo- 
scalar vertex functions in the momentum space in the Landau gauge and applied the RI/MOM 
scheme [^. In the fit of the vertex functions we explicitly use the low-mode contribution to the 
chiral condensate 

with eigenvalues A,- and absorb its valence quark mass dependence. By converting zfl^^'^^ in the 
massless limit into the value of MS perturbatively, we obtain a preliminary value Z^^(2GeV) = 
0.815(8). 

On our lattice with spatial extent L, the lightest pion gives mj^L w 2.8. Finite size effect (FSE) 
may therefore be sizable. We estimate this effect using the analytic result [||] from a combination 
of the resummed Luscher's formulae and an one-loop ChPT. There is an additional finite size due 
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Figure 1: Effective pion mass plots for the four lightest quark mass. Open symbols indicate signals from 
conventional treatment while filled symbols are from low-mode-averaging. 



to the fixed topology in our numerical simulation [[lO|]. Based on the discussion in [11], we make 
a correction using the result of one-loop ChPT and the numerical data of topological susceptibil- 
ity [|I2|] determined on the same lattice configurations. 



3. Chiral extrapolation {Nf = 2) 

We address the validity of the NLO ChPT prediction which is commonly used in the chiral 
extrapolation of spectrum data on the lattice In the framework of Nf = 2, pion mass and decay 
constants are expanded in terms of x = ABniq/ (47r/)^ as 

m^/niq = 2B{\ + jx\nx)+cj,x, (3.1) 
fn = /(I — xlnx) + C4X (3.2) 

to NLO {i.e. one-loop level or i^{x)), where B and / are the tree level low energy constants (LECs), 
and C3 and C4 are related to the one-loop level LECs li, and 1^^. At NLO, these expressions are 
unchanged when one replaces the expansion parameter x by x = 2m^/{4nf)^ or ^ = 2m\/{AKf-jtY, 
where m\ and fj^ denote those at a finite quark mass. Therefore, in a small enough pion mass region 
the three expansion parameters should describe the lattice data equally well. 

Three fit curves (x-fit, x-fit and t, -fit) for the three lightest pion mass points {rriji < 450 MeV) 
are shown in Figure § as a function of For all fits, the horizontal axis is appropriately rescaled 
to give m| using the obtained fit curves. 

From the plot we observe that the different expansion parameters describe the three lightest 
points equally well; the values of ;^^/dof are 0.30, 0.33 and 0.66 for x-, x- and <^-fits, respectively. 
In each fit, the correlation between m|/m^ and fj^ for common sea quark mass is taken into account. 
Between the x- and x-fit, all of the resulting fit parameters are consistent. Among them, B and / are 
also consistent with the i§ -fit. This indicates that the NLO formulae successfully describes the data. 
In Figure |[, results of B (upper panel) and / (bottom panel) for different fits are plotted for different 
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Figure 2: Comparison of the chiral fits including the NLO terms for m\/m^ (left) and /^^ (right) obtained 
in the Nf — 2 calculation. Fit curves for the three lightest data points obtained with different choices of the 
expansion parameter (jc, x and ^) are shown as a function of m^. 
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Figure 3: Results of fit parameters B (top) and / (bottom) as functions of the upper limit of the fit range. In 
each panel, circle, square and diamonds are obtained with fit parameters x, x and ^ . Results with ;i;^/dof < 2 
are plotted. 



pion mass points. As seen in the figure, the agreement among the different expansion prescriptions 
is lost when we extend the fit range to include the 4th lightest data point which corresponds to niji ~ 
520 MeV. We, therefore, conclude that for these quantities the NLO ChPT may be safely applied 
only below « 450 MeV. 

Another important observation from Figure ^ is that only the ^ -fit reasonably describes the 
data beyond the fitted region. With the x- and x-fits the curvature due to the chiral logarithm is too 
strong to accommodate the heavier data points. In fact, values of the LECs with the x- and x-fits are 
more sensitive to the fit range than the ^ -fit. This is because /, which is significantly smaller than 
fjz of our data, enters in the definition of the expansion parameter. Qualitatively, by replacing niq 
and / by m| and the higher loop effects in ChPT are effectively resummed and the convergence 
of the chiral expansion is improved. 
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Figure 4: Comparison of the Nf = 2 results from the NLO fit and the NNLO fit with ^ . Black pluses denote 
reference points from phenomenological estimations. 



We then extend the analysis to include the NNLO terms [|T^]. Since we found that only the E,-f\t 

reasonably describes the data beyond m-^ ~ 450 MeV, we perform the NNLO analysis using the t, - 

expansion. Although we input phenomenological estimate of LECs l\ and Ij, we find our fit result 

is insensitive to their uncertainties. We extract the LECs of ChPT, i.e. the decay constant in the 

chiral hmit /, chiral condensate £ = Bf/2, and the NLO LECs /f''^' = -c^/B + ln{2V2Kf /m,,+ f 

and if'^^ = €4// + ln{2\/2nf /m-^+y. For each quantity, a comparison of the results between the 

NLO and the NNLO fits is shown in Figure 0. In each panel, the results with 5 and 6 lightest data 

points are plotted for the NNLO fit. The correlated fits give ;^^/dof = 1.94 and 1.40, respectively. 

For the NLO fits, we plot results obtained with 4, 5 and 6 points to show the stability of the fit. The 

X^/dof is less than 1.94. The results for these physical quantities are consistent within either the 

NLO or the NNLO fit. On the other hand, as seen for if^^^ most prominently, there is a significant 

disagreement between NLO and NNLO. This is due to the large NNLO contributions to the terms 

which are proportional to C3 and C4, respectively. 

We quote our final results for the Nf = 2 calculation from the NNLO fit with all data points: 
/=111.7(3.5)(1.0)C^;»)MeV,r^^(2GeV) = [235.7(5.0)(2.0)Ci2.7)^gyp^;-p^^^ 

and If^' = 4.12(35)(30)(+^|J), where m+ = 139.6 MeV. From the value at the neutral pion mass 
nijfi = 135.0 MeV, we obtain the average up and down quark mass m„d and the pion decay con- 
stant as m^{2 GeV) = 4.452(81)(38)(+227) MeV and fj^ = 119.6(3.0)(1.0)(+^:^) MeV. In these 
results, the first error is statistical, where the error of the renormalization constant is included in 
quadrature for Z'/^ and m„d. The second error is systematic due to the truncation of the higher 
order corrections. For quantities carrying mass dimensions, the third error is from the ambiguity in 
the determination of tq. We estimate these errors from the difference of the results with our input 
ro = 0.49 fm and that with 0.465 fm. The third errors for if^^^ and if^^'^ reflect an ambiguity of 
choosing the renormalization scale of ChPT {47if or 47ifji). 
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4. Chiral extrapolation {Nf = 2+1) 

Since we found in the two-flavor calculation that the NNLO ChPT formulae can nicely fit our 
data even in the kaon mass region if one uses the -expansion, we apply the same strategy for our 
2 + 1 -flavor analysis. 

As functions of E,:,^ = m^/{4Kfji)^ and ^ri, = fn\J {^T^fn)^, where r\s is the unphysical strange- 
strange meson, chiral expansions are expressed as 

ml/niud = 2Bo(l +M\^,, ^^-U)) + a^^l + a^^^^r,. + a^^r,^ (4.1) 

m\lma, = 2Boil+M''{^^,^^-L';)) + afU^n-U) + (^U(U-^^) (4-2) 
A = /o(l+F'^(^.,^,.,;Ln) + i8f<§2 + ^2^^^^^^ + ^3-^2 ^43) 

fK = /o(l+F^(^.,^r,.;LD)+i8f§.(^.-^,J+i8f^r,,,(^,.-^.), (4.4) 

where rriav = jims + rriud) and a^^^'^ and jS^^'^ are NNLO unknown parameters. Functions M^, M^, 
and contain NLO contributions and loop contributions at NNLO, whose expressions are too 
involved to present here [p^]. Among relevant SU(3) LECs L^-L^, we use values L[ = (0.38 ± 
0.18) • 10-^ L$ = (1.59 ±0.15) • 10-^ L5 = (-2.91 ±0.32) • lO^^ andL^ = (-0.49 ±0.24) • 10"^ 



(defined at /i = 770 MeV) from a phenomenological estimate [15] and determine others by a fit. 
Thus, the chiral extrapolation with ( |4. l[ )-(|44l) contains 16 fit parameters in total. We fit 
m^/niav, fn and fx simultaneously taking the correlation between these quantities at the same sea 
quark mass {mud,ms) into account. By using all data points, x^/dof = 1.42 is obtained. Data after 
the finite size corrections are shown in Figure |5| as a function of Different symbols correspond 
to the pion data (m^/mud and fj^) and the kaon data {rriK/mav and fx) while the filled (open) pattern 
represent a fixed lighter (heavier) strange quark mass, which is accompanied by the solid (dashed) 
curves. 

Extrapolating the data to the physical point rrij^ = 135.0 MeV, rriK = 495.0 MeV and fj^ = 130.7 
MeV, we obtain prehminary results m^(2 GeV) = 3.799(68) MeV, mp(2 GeV) = 1 14.6(2.0) MeV, 
fn = 121.5(4.1) MeV, fx = 148.3(4.7) MeV undfK/fn = 1.220(10), where the errors are statisti- 
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cal only. Further studies to determine LECs and to analyze their flavor dependence are planned as 
well as the estimation of systematic errors. 

Numerical simulations are performed on Hitachi SRllOOO and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 07-16). This work is supported in part by the Grant-in-Aid of the Ministry 
of Education (Nos. 18340075, 18740167, 19540286, 19740121, 19740160, 20025010, 20039005, 
20340047, 20740156), the National Science Council of Taiwan (Nos. NSC96-2112-M-002-020- 
MY3,NSC96-2112-M-001-017-MY3, NSC97-2119-M-002-001) and NTU-CQSE (No. 97R0066- 
69). 
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